Heat conduction: a telegraph-type model with 
self-similar behavior of solutions 
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For heat flux q and temperature T we introduce a modified Fourier-Cattaneo law q t + lj — —kT x . 
The consequence of it is a non-autonomous telegraph-type equation. This model already has a 
typical self-similar solution which may be written as product of two travelling waves modulo a 
time-dependent factor and might play a role of intermediate asymptotics. 
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It is well known, that the heat equation propagates perturbations with infinite velocity. For this contradiction a 
possible answer is the telegraph equation which is " obviously hyperbolic" . In general, one usually forgets about the 
other fundamental properties of parabolic heat equations: the existence of self-similar solutions (e.g., the Gaussian 
kernel or fundamental solution) and the attracting nature of these special solutions (intermediate asymptotics). 

Is is easy to show, that the telegraph equation -see (5) -has no self-similar solutions e.g., solutions of the form 
t~ a f(x/t@) and that even asymptotic self-similarity property is lacking: no solutions of form g(t) ■ f(x/w(t)) with 
43 I g ~ t- a and h ~ t p for t » 1. 
©H" Since the telegraph equation (possibly with reaction terms) supposed to be relevant not only in heat conduction 
but also in various diffusion processes, the lack in self-similarity might be a bad sign for the adequacy of the model. 
In addition, in diffusion and heat theory various physical quantities, like the fluxes, have to be continuous; therefore 
the solutions of this equation cannot be " too bad" . 

According to Gurtin and Pipkin [lH3|, the most general form of the flux in linear heat conduction and diffusion 
related to the flux q expressed in one space dimension via an integral over the history of the temperature gradient 



on ' 

ON ' where Q(t — t') is a positive, decreasing relaxation function that tends to zero as t — t' — > oo and T{x,t) is the 

\ temperature distribution. 

. There are two notable relaxation kernel functions: if Q(s) — kS(s) where 5(s) is a Dirac delta "function", then 

§ : , = - fe ^M (2) 
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is the Fourier law. 

If we define the kernel as Q(s) = ^e~ s ' T where s =t — t and k is the constants of the effective thermal conductivity 
5_i ' we get back to the well-known Cattaneo [4J heat conduction law which reads: 

dq dT(x,t) 
T d-t +q= - k ^x— (3) 

fc2 p —s/t 



There is a third kind of relaxation kernel, the " Jeffrey- type" , which applies both Q(s) — kiS(s) + 
The energy conservation law 

where 7 is the heat capacity, gives the heat equation with (1) and the telegraph equation with (2): 

d 2 T(x,t) 1 dT(x,t) _ 2 d 2 T(x,t) 

Ot 2 + t 8t ~° 8x 2 ' [ > 



where c = yfk/r^ is the propagation velocity of the transmitted heat wave. The flux q satisfies the same equation. 

The thermal diffusivity k = fc/7 can be defined as the ratio of the effective thermal conductivity k and the heat 
capacity 7. This equation describes a physical process with a well defined relaxation time r. 



2 



The telegraph equation ([5]) can be derived in various transport systems. First it was derived by Kirchhoff to 
describe the voltage on a realistic electrical transmission line with distance and time, this derivation can be found 
in any textbook on electrodynamics. Later Goldstein [j| derived the telegraph equation from a generalized random 
walk for diffusion. Okubo [6] summarized the knowledge of the telegraph equation till 1971 and presented a new 
ad-hoc derivation for turbulent mixing his starting point being the Navier-Stokes system. In all these derivations the 
t relaxation time had a well-established physical meaning presenting the time-scale of the physical process. Joseph 
and Preziosi collected basically all the relevant works done on heat waves connected with the telegraph equation till 
1990 [2J, |3|]. On the other side the telegraph equation was intensively investigated from the mathematical point of 
view as well Q ■ Different kind of properties of these models can be analyzed as in @ . 

In the present communication we introduce a new kernel which somehow interpolates the Dirac delta and the expo- 
nential kernel having the main properties of both. The Q(s) — 1/s is such a function: it is singular at the origin and 
has a short range of decay for I > 1. Let's consider the following relaxation kernel 

where k is the effective thermal conductivity, r is a relaxation time and I > 1 is a parameter, — t' + uj is just a time 
shift which is needed to regularize the expression. 



Using the general form of heat flux ([T]) we get 



kTl dT ^-dt'. (7) 



(t-f + uj) 1 dx 



dq = k fT_\ l dT(x,t) 



One has 

dt \ujJ dx 

1 [* 1 kr l dT(x,t') 

l J_ go t-*+ u {t-i> + uy dx l8J 

A formal application of integral mean theorem to the second term on the right hand side and the definition of q leads 
to a new phenomenological law 



dq , (T\i dT(x,t) I 

The additional energy conservation law is still ((4]) and from this and the last equation we obtain 

7 fuj\' d 2 T(x,t) 



(-) 

k\rJ dt 2 

7 fw^l I dT(x,t) _ d 2 T(x,t) 
t-t" + uj dt ~ dx 2 ' 



(10) 



For a better transparency let's call e = ¥ (— ) arid a = ¥ (— ) • I The physical meaning of e is still the thermal 
diffusivity multiplied by a scaling constant which is the renormalized relaxation time (the ratio of an ordinary time 
shift uj and a well defined relaxation time r). The exponential I is a real number which describes the non- locality 
in time which we may call memory effects of the heat conduction phenomena. Larger I means shorter memory. 
The physical meaning of a is approximately the thermal diffusivity multiplied by another time-scaling factor. In the 
following we will see that the role of a, e or I will be crucial in the structure of the solutions. At last we introduce a 
new time variable t — t — t" + uj. Now our telegraph- type equation reads 

d 2 T{x,t) adT{x,t) _ d 2 T(x,t) 

e ~d£ + 1 dt ~ dx 2 ' 

Note, that the a/t factor appearing in front of the first time derivate makes the equation time- reversible, which cannot 
be true for diffusion or heat propagation processes, at the same time the a/t factor makes the equation irregular at 
the origin. To avoid these problems, we may shift the pole to a negative time value - in practical calculations using 
the a/(t + t) where r still can be any kind of relaxation time with well-founded physical interpretation. Physically 
it is clear, if a process has a well-defined time-scale than the reverse process cannot run back in time more than the 
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physically relevant time. Now, in this sense, for positive time (t > 0) we may use the equation to describe diffusion-like 
processes. Our deeper investigation clearly shows, that no other exponent of t as one can have self-similar solution. 

If we consider as a non-linear wave equation we may investigate wave properties like dispersion phenomena. In- 
serting the standard plain wave approximation T(x, t) — e l i kx + UJt ) into the dispersion relation and the attenuation 
distance can be obtained. These are the followings: 

-1/2 

Cj 12 I I flY \ _ 1 2t 1 




, ~ ......... , ■ (12) 

Re(k) Ve I V W ) Im(k) d v p 

Our telegraph-type equation is time dependent, hence both the dispersion relation and the attenuation distance have 
time dependence. Note, that v p has a very weak time-dependence, basically only till t < I. As an other interesting 
point is that the phase velocity does not depend on the angular velocity which is the same as for the ideal wave 
equation, so our equation has no dispersion. So in this sense our equation is very similar to the wave-equation, which 
is hyperbolic. The properties of the attenuation distance is even more interesting, it is divergent in time and has a 
1/uj angular frequency. However if we let the angular frequency and the time to go infinite with the same speed than 
the attenuation distance has a strong decay. Which is like the skin-effect when high frequency electrons can only 
propagate on the surface of a metal. 

We are looking for solution of (fTTj) of the form 

T(x,t) = *""/( :=t- a f( V ). (13) 

The similarity exponents a and /3 are of primary physical importance since a represents the rate of decay of the 
magnitude T(x,t) 1 while /3 is the rate of spread (or contraction if j3 < ) of the space distribution as time goes on. 
Substituting this into (TTT]) we have 

f(ri)^- a - 2 [eaP - e(3(-a - (3 - 1) - (3a\ + 
f( V )t- a - 2 [-ea(-a - 1) - aa] = /"fa)*-"" 2 * (14) 

where prime denotes differentiation with respect to r\. 

One can see that this is an ordinary differential equation(ODE) if and only if a + 2 = a + 2/3 (the universality 
relation). So it has to be 

(3 = 1 (15) 

while a can be any number. The corresponding ODE we shall deal with is 

f'\ri)[eri 2 - 1] + f'(rj)r](2ea + 2e - a) + f(n)a(ea + e - a) = 0. (16) 

In pure heat conduction-diffusion processes-no sources or sinks-the heat mass is conserved: the integral of iXcc, t) 
with respect to x does not depend on time t. For T(x, t) this means 

J T(x, t)dx = t~ a J f(j)dx = t- a+1 J f{n)dn = const (17) 

if and only if a = 1. We are going to investigate this case only. Plainly (flT)]) can be written as 

[efv 2 /)" = a( V f)' (18) 
which after integration and supposing f(r)o) — for some 7/0 gives 

df andrj 



f m 
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(19) 



From this equation we can obtain two qualitatively different solutions. The one which is globally bounded and positive 
in the domain {(x, t) : 1 — erj 2 > 0} and has the form 

f = (l-e V 2 )f- 1 (20) 
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FIG. 1: Eq. (21) thick solid line is for I — 6.2 and the thin dashed line is for I = 4.1. 




FIG. 2: The solution (|21j) for the parameter I 



where (/)+ = max(/,0). See Fig. ©. 
The corresponding self-similar solution is 



1 / x 2 



(21) 



This solution is positive in the cone t 2 > ex 2 and is zero outside of it, see Fig. ([2]). Note, that only the x > and 
t > quarter of the plane is presented, because of it has physical relevance. 

On the (x,t) plane there are two fronts x(t) = ±-Jg separating these domains. Because the function T(x,t) not 
always has continuous derivatives entering to (1111) we have to make clear what we mean under " solution" . Having 
in mind the physical background, we ask the continuity of T t ,T x ,q t and q x so that in (2) and (3) all functions were 
continuous. In our case this means that 



a I — 

2e~ ~ ~ 2 



(22) 



has to be greater than 1, i.e. a/e — I > 4, which we shall suppose further on. If the second derivatives are not 
continuous (4 < I < 6) we understand the solution in the sense of distributions. If I > 6 the solution is classical. On 
Fig. 1. we compare the solutions with I = 4.1 and / = 6.2. the thick solid line represents the solution for I — 6.2 and 
the thin dashed line however shows the solution for I = 4.1. 
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Remark 1. The solution (121) is of source-type i.e. lim^o T(x, t) — K5(x), where <5 is the Dirac measure, K > 0. 
One can calculate the second initial condition lim^o T t (x, t) too. 

Remark 2. One can write (|21|) in the form of product of two traveling waves propagating in opposite direction 
(divided by a time- factor): 

T(x,t) = -? ZI (t-s/ix)i-\t + s/ix)l- 1 (23) 

which is a new-type of purely hyperbolic wave; the typical solution of the wave equation is the sum of two such waves: 
g(x — ct) + g(x + ct). 

It is known that an another possible answer to contradiction connected with the infinite speed of propagation is 
the nonlinear Fourier law (r = 0, k = k^T" 1 ^ 1 in (2)) which leads to a nonlinear heat equation 

T t = (T m ) xx , m > 1. (24) 

In [lol ] Zeldovich and Kompaneets have found the fundamental solution T\ of this equation which we write in the 
following form: 

T rn-1 = t _a(m-l) (^2 _ £2^-2^ = 



-{AtP -Bx)+(AtP +Bx)+, (25) 



where A is constant and 



a = j3 = —!—, B 2 = - m \ . (26) 
m + 1 ' 2m(rn + 1) v ; 

One can see that this solution has bounded support in x for any t > which is a hyperbolic property. Using comparison 
principle for such equations one can show this finite speed property for any initial condition having compact support. 
However, the fronts are not straight lines: x(t) = i-gf 9 , /3 < 1 so the speed of propagation x(t) goes to zero if t goes 
to infinity. One can also see that T\ is of source-type: Ti(x, 0) = KiS(x). 

The most intrinsic property of T\ is that it plays the role of intermediate asymptotic: any solution of (|24p corre- 
sponding to the initial datum t(x,0) with J t(x j 0)dx — K\ converges to T\ as t — > oo. This was conjectured earlier 
but was shown only in 1973 by Sh. Kamin, see 

It would be important and interesting to understand whether or not our special solution T(x,t) had this attractor 
property. If "yes", in what sense: we recall that there is a second initial condition too. 

In summary. - We introduced a new phcnomcnological law for heat flux which in some sense " interpolates" between 
Fourier and Cattaneo laws. The consequence of it is a non-autonomous model, a telegraph- type partial differential 
equation. It already has, unlike the classical telegraph equation, self-similar solutions, the presence of which is desirable 
in the theory of heat propagation free from sources and absorbers. 

One of us (R.K.) would like to thank Prof. P. Rosenau for his stimulating discussion. 
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